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Introduction 


While  running  has  been  shown  to  reduce  disease  risks  and  promote  a  generally  healthy  lifestyle 
in  uninjured  people,  very  little  running-specific  r  esearch  is  available  pertaining  to  the  amputee 
population.  The  little  existing  a  mputee  running  literat  ure  primarily  involves  running  with 
prostheses  designed  for  every  day  wear,  which  are  typically  prescribed  and  aligned  to  perform 
optimally  during  standing  and  walking.  Further,  these  studies  have  used  biomechanical  models 
designed  for  the  intact  limb  to  estimate  joint  kinetics  (f  orces  and  moments)  usin  g  an  inverse 
dynamics  approach.  This  approach  estimates  distal  joint  kinetics  and  uses  these  calculations  to 
estimate  more  proxima  I  joint  kinetics.  C  onsequently,  inverse  dynamics  est  imations  rely  on 
accurate  estimation  of  the  ankle  joint  as  erro  rs  will  be  p  ropagated  and  inflated  with  more 
proximal  calculations.  The  previous  stud  ies  on  amputee  running  have  not  validated  the 
methodology  used  for  joint  kinet  ic  measurements  with  running-specific  prostheses,  which  can 
potentially  prove  to  be  erroneou  s.1'3  These  limitations  call  for  systematic  research  an  d 
development  of  validated  models  for  running-sp  ecific  prostheses.  This  will  lead  to  improved 
prosthetic  designs  that  will  allow  clinicians  to  provide  evidence-based  exercise  prescriptions  to 
amputees,  enabling  th  em  to  comfortably  and  efficiently  r  un.  The  objective  of  the  propose  d 
study  is  to  develop  and  validate  a  model  with  optimal  set-  up  of  reflective  markers  used  in  3D 
gait  analysis,  producin  g  minimal  errors  in  in  verse  dynamics  calculations.  The  long-term 
objective  of  this  project  is  to  understand  the  biomechanical  and  physiological  consequences  of 
amputation,  to  develop  an  optimal  design  of  activity-independent  lower-extremity  prosthesis, 
and  to  help  clinicians  prescribing  appropriate  prosthesis  and  exercise  regimes  to  people  with  a 
lower  extremity  amputation. 

Please  see  Appendix  II:  Manuscript  Draft  1  for  a  more  detailed  introduction. 

The  current  project  was  approved  for  funding  over  18  months  beginning  1  August,  2009  and  to 
be  completed  by28F  ebruary,  2011.  A  no-cost  extensio  n  was  approved  due  to  issue  s  in 
prosthesis  procurement  caused  by  the  prosthetic  manufacturing  companies.  The  extension  was 
approved  through  28  February,  2012. 

The  purpose  of  this  do  cument  is  to  provide  details  of  the  study  to  satisfy  the  Final  Report 
requirement. 
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Body 

The  approved  Statement  of  Work  proposed  the  following  timeline  (Table  1): 


Table  1.  Timeline  for  approved  project. 


Months  1-6 

Months  7-12 

Months  13-18 

Specific  Aim  #1 :  Development  and  validation  of  a  model  with  unique  optimal  marker 
placements  for  specific  running  prosthesis  designs 

Formulate  program  for  data  analysis 

X 

MTS  testing  of  running  specific 
prostheses  (12  prostheses) 

X 

Validation  of  model 

X 

Analysis  of  MTS  data  to  determine  model 

X 

Specific  Aim  #2:  Determination  of  the  resu 
running  prosthetic  designs 

tant  optimal  marker  placement 

for  all  tested 

Determine  model  with  optimal  marker 
placement  for  across  designs 

X 

Methodology 

A  biomechanical  model  was  developed  usin  g  motion  a  nalysis  of  running-spe  cific 
prostheses  in  a  material  testing  system  (MTS,  Eden  Prairie,  MN).  Four  running-specific 
prosthesis  designs  (Figure  1)  were  tested  for  this  project  including  the  1 E90  Sprinter  (OttoBock 
Inc.),  Flex-Run  (Ossur),  Cheetah  ®  (Ossur)  and  Nitro  Running  Foot  (Freedom  Innovations). 
These  prostheses  were  chosen  because  they  are  the  m  ost  commonly  prescrib  ed  running- 
specific  prostheses  currently  available  on  the  market.  Three  different  stiffness  categories  were 
also  tested  for  each  prosthetic  de  sign  to  identify  whether  prosthetic  stiffness  affe  cts  optimal 
marker  placement.  Stiffness  categories  were  chosen  to  reflect  a  common  range  of  stiffnesses 
that  might  be  prescribed. 

Each  prosthesis  was  placed  in  the  MTS  between  two  load  cells  (Bertec  PY6,  Columbus, 
OH)  in  a  neutral  alig  nment.  Neutral  align  ment  was  defined  according  to  t  he  specific 
manufacturers’  recommendations  for  prosthesis  alignment.  The  load  cells  captu  red  data  at 
1 ,000  Hz.  The  prosthe  ses  were  cyclically  loaded  for  ten  cycles  with  axial  forces  u  p  to  2,500  N 
to  simulate  peak  vertical  forces  commonly  obs  erved  during  running  (approximately  t  hree  times 
the  body  weight  of  a  75  kg  person).  The  load  cells  measured  the  force  and  moment  at  the  point 
of  load  application  proximal  to  the  prostheses  (applied  load)  and  the  reaction  forces  distal  to  the 
prostheses  (ground  reaction  forces). 

Reflective  markers  were  placed  at  2  cm  intervals  along  the  lateral  aspect  of  the  keel  of 
each  running-specific  prosthesis  (Figure  2).  Reflective  markers  were  also  placed  o  rthogonally 
on  the  anterior,  lateral,  and  media  I  aspect  of  the  “head”  of  the  prosthesis,  at  th  e  point  of 
connection  to  the  socket  or  pylon,  in  order  to  define  the  local  coo  rdinate  system  of  the 
prosthesis.  Three  additional  markers  were  placed  along  the  midline  of  each  prosthesis  to  define 
a  plane  to  which  the  keel  markers  were  projected  for  further  analysis.  An  8-camera  motion 
capture  system  (Vicon,  Oxford,  UK)  with  a  capt  ure  frequency  of  500  Hz  was  used  to  collect  the 
3-D  positional  data  of  the  markers  during  each  trial.  Two  consecutive  projected  center  lin  e 
markers  defined  individual  segments  of  the  prosthesis  (a  ssumed  to  be  rigid)  and  consecut  ive 
segments  shared  a  co  mmon  marker.  The  joint  between  these  segments  was  assumed  as  a 
hinge  joint.  Standard  inverse  dynamics  calculations  (equations  1  and  2)  were  used  to  estimate 
the  force  and  torque  transfer  from  the  ground  reaction  force,  through  the  defined  prosthesis 
segments,  and  to  the  point  of  load  application  proximal  to  the  prosthesis. 
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a.  Freedom  Innovations 
Nitro 


b.  Ossur  Flex-Run 


c.  Ossur  Cheetah 


d.  Ottobock  Sprinter 
1E90 


Figure  1.  Prostheses  used  for  mechanical  testing. 


+  [1] 
iMf  +■  ^  [2] 

where  Fi  and  Fi+1  are  forces  acting  on  link  /'  at  joints  /  (distal)  and  i+1  (proximal),  respectively;  Mi, 
Mi+1  are  moments  exerted  on  link  /  at  joints  /  and  i+1,  respectively;  n  and  ri+i  are  radii  from  the 
COM  of  link  /  to  the  joint  centers  /  and  i+1,  respectively;  g  is  acceleration  due  to  gravity;  [/,]  is  the 
matrix  of  inertia;  and  and  i?,  are  vectors  of  the  angular  velocity  and  acceleration  for  link  /, 
respectively. 


Please  see  Appendix  II:  Manuscript  Draft  1  for 
additional  details  and  figures  related  to  the  study 
methodology. 

Research  Accomplishments 

Specific  Aim  #1 


MTS  testing  of  running  specific  prostheses:  Completed 
The  experimental  setup  was  finalized  and  MT  S  testing 
was  performed  for  t  he  existing  running-specific 
prostheses.  A  representative  example  of  a  prosthesis 
set  up  in  the  MTS  machine  is  shown  in  Figure  2 
Reflective  markers  were  placed  a  long  the  ke  el  of  the 
prosthesis  in  2  cm  increments  and  force  transducers 
are  present  at  the  base  and  top  of  the  experimental 
setup  in  order  to  measure  forces  an  d  moments  at  the 
input  (top)  and  “ground”  level. 


Figure  2.  Running-specific 
prosthesis  (Ossur  Flex-Run)  setup 
in  the  MTS  machine. 


The  difference  between  force  and  moment  values  at  the  point  of  load  application  from 
the  estimated  inverse  dynamics  calculations  and  the  directly  measured  values  from  the  top  load 
cell  was  considered  model  error.  Force  and  mo  ment  estimations  were  made  with  every 
combination  of  remaining  markers  giving  a  resultant  error  value  for  ea  ch  combination.  Errors 
were  calculated  for  ea  ch  loading  cycle  as  ro  ot  mean 
squared  error  (RMSE)  and  normalized  RMSE 
(NRMSE),  respectively. 


Formulating  program  for  data  analysis:  Completed 

Once  raw  experimenta  I  data  were  obtained  from  the  MTS  testing,  we  were  able  to  beg  in 
formulating  the  data  analysis  progr  am  in  the  Matlab  programming  la  nguage  and  validate  the 
proposed  model.  Because  the  project  proposed  a  new  model  and  met  hod  to  analyze  running- 
specific  prostheses,  ea  ch  stage  of  the  progra  m  development  required  validation  to  ensure 
proper  measurements  and  calculations.  Consequently,  this  was  a  lengthy  progress  involving  a 
large  amount  of  troubleshooting.  The  data  analysis  program  was  completed  so  the  results  could 
be  used  to  validate  the  marker  model.  Please  see  the  data  generated  from  this  programming  in 
the  Reportable  Outcomes  section  of  this  document. 

Validation  of  model:  Completed 

The  programming  and  model  validation  were  completed  for  use  in  the  analysis  of  the  MTS  data 
to  determine  the  final  marker  model. 

Analysis  of  MTS  data  to  determine  model:  Completed 

Determining  the  final  marker  model  for  each  s  pecific  prosthesis  design  was  completed. 
Completion  of  this  task  was  delayed  beyond  the  originally  proposed  timeline  due  to  procurement 
issues  with  the  prosthetic  companies.  These  issues  were  resolved  and  the  task  was  completed. 
This  task  completes  the  goals  set  for  Specific  Aim  1 . 

Specific  Aim  #2 

Determine  model  with  optimal  marker  placement  for  all  prosthesis  designs:  Completed 
This  task  was  proposed  for  Months  13-18  of  the  project  and  required  the  completion  of  tasks  in 
Specific  Aim  1 .  This  task  was  delayed  due  to  prosthesis  procurement  issues.  However,  th  e 
task  was  completed  as  proposed. 


Results 


Data  are  presented  for  anteroposterior  (AP)  force,  vertical  force,  an  d  flexion  momen  t 
values  throughout  the  cyclical  load  ing  task.  Fi  gure  3  compares  the  directly  measured  values 
with  those  estimated  from  thousands  of  d  ifferent  combinations  of  marker  placements  on  t  he 
prosthesis.  Figure  4  displays  the  r  aw  error  between  the  estimated  force  and  moment  values 
and  the  directly  measu  red  values  for  each  of  the  marker  placement  combinations.  Figure  5 
shows  the  NRMSE,  representing  the  difference  (in  percent)  between  the  directly  measured  (via 
load  cell)  a  nd  estimated  (via  inverse  dynamics)  proximal  forces  and  moments  calculated  for 
each  marker  combinations. 

The  Freedom  Innovations  Nitro  prosthesis  had  a  maximal  RMSE  range  of  0.26  N  (AP 
force),  4.45  N  (vertical  force),  and  1.02  Nm  (fie  xion  moment)  and  a  maximal  NRMSE  range  of 
0.02%,  0.17%,  and  0.86%  for  AP  fo  rce,  vertical  force,  and  flexion  moment,  respectively  across 
all  stiffness  categories  and  all  tested  combinations  of  markers.  The  Ossur  Flex-Run  prosthesis 
had  a  maximal  RMSE  range  of  4.37  N  (AP  force),  5.88  N  (vertical  force),  and  1 .05  Nm  (flexion 
moment)  and  a  maximal  NRMSE  range  of  0.37%,  0.28%,  and  0.56%  for  AP  force,  vertical  force, 
and  flexion  moment,  respectively  across  all  st  iffness  categories  and  all  tested  combinations  of 
markers.  The  Ossur  Cheetah  prosthesis  had  a  maximal  RMSE  range  of  0.99  N  (AP  force),  9.38 
N  (vertical  force),  and  0.73  Nm  (flexion  moment)  and  a  maximal  NRMSE  range  of  0.12%, 
0.44%,  and  0.53%  for  AP  force,  vertical  for  ce,  and  flexio  n  moment,  respectively  across  a  II 
stiffness  categories  and  all  tested  combinations  of  markers.  The  Ottobock  1E90  prosthesis  had 
a  maximal  RMSE  rang  e  of  0.48  N  (AP  force),  7.54  N  (ve  rtical  force),  and  0.54  Nm  (flexion 
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moment)  and  a  maximal  NRMSE  range  of  0.07%,  0.35%,  and  0.31%  for  AP  force,  vertical  force, 
and  flexion  moment,  respectively  across  all  st  iffness  categories  and  all  tested  combinations  of 
markers. 

Figures  6-8  show  the  average  RMSE  values  for  AP  force,  vertical  force,  and  flexion 
moment,  respectively,  for  each  prosthesis  across  the  number  of  markers  on  the  prosthesis.  All 
tested  combinations  with  the  number  of  markers  indicated  were  averaged  to  generate  each  data 
point. 

Cumulatively,  these  da  ta  indicate  little  differen  ce  in  kinetic  calculat  ions  between  the 
directly  measured  values  and  any  marker  placement  or  combination  of  markers.  Consequently, 
a  single  marker  can  be  placed  on  any  point  of  a  running-specific  prosthesis  to  achieve  accurate 
proximal  kinetic  estimations. 

Please  see  Appendix  II:  Manuscript  Draft  1  for  a  more  detailed  report  of  the  study  results. 
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Flexion  Moment 


150 
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Figure  3.  AP  force,  vertical  force,  and  flexion  moment  curves  for  cyclical  loading.  Thick  blue 
lines  represent  the  directly  measured  values  fr  om  the  upper  load  cell.  Thin  lines  represent 
calculated  values  from  each  different  combination  of  mark  ers.  Exemplar  data  is  from  th  e 
Ossur  Flex  Run  categ  ory  3  prosthesis.  Other  tested  prostheses  a  nd  stiffness  categories 
showed  similar  results. 


AP  Force 


Vertical  Force 


Flexion  Moment 


^0  50  100 

%  Loading  Cycle 


Figure  4.  Average  AP  force,  vertical  force,  a  nd  flexion  moment  error  curves  for  t  he  loading 
cycle.  Each  curve  repr  esents  the  difference  between  the  directly  measured  values  from  the 
upper  load  cell  and  calculated  values  from  each  combination  of  marke  rs  on  the  prosthesis. 
Exemplar  data  is  from  the  Ossur  Flex  Run  category  3  prosthesis.  Other  tested  prostheses  and 
stiffness  categories  showed  similar  results. 
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AP  Force  Vertical  Force  Flexion  Moment 


#  of  Markers  on  Prosthesis  #  of  Markers  on  Prosthesis  #  of  Markers  on  Prosthesis 

Figure  5.  Normalized  root  mean  square  error  (NRMSE)  for  each  combination  of  markers  for 
AP  force,  vertical  for  ce,  and  flexio  n  moments  throughout  the  load  ing  cycle.  Each  dot 
represents  the  NRMSE  value  for  a  particular  combination  of  markers.  The  x-axis  shows  the 
number  of  markers  on  the  prosthesis  for  the  particular  combination.  Exemplar  data  is  from  the 
Ossur  Flex  Run  categ  ory  3  prosthesis.  Other  tested  prostheses  a  nd  stiffness  categories 
showed  similar  results. 
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Figure  6.  Average  anteroposterior  (AP)  force  root  mean  squared  error  (RMSE) 
for  each  prosthesis  across  the  number  of  mark  ers  on  the  prosthesis.  All  tested 
combinations  with  the  number  of  markers  indicated  were  averaged  to  generate 
each  data  point.  Error  bars  represent  ±1  standard  deviation  of  all  marker 
combinations  tested  for  the  number  of  markers  shown. 
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Figure  7.  Average  vertical  force  root  mean  squared  err  or  (RMSE)  for  each 
prosthesis  across  the  number  of  markers  on  the  pro  sthesis.  All  tested 
combinations  with  the  number  of  markers  indicated  were  averaged  to  generate 
each  data  point.  Error  bars  represent  ±1  standard  deviation  of  all  marker 
combinations  tested  for  the  number  of  markers  shown. 
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Figure  8.  Average  flexion  mome  nt  root  mean  squared  error  (RMSE)  for  each  prosthesis 
across  the  number  of  markers  on  the  prosthesis.  All  tested  combinations  with  the  number  of 
markers  indicated  were  averaged  to  generate  each  data  point.  Error  bars  repr  esent  ±1 
standard  deviation  of  all  marker  combinations  tested  for  the  number  of  markers  shown. 
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Key  Research  Accomplishments 


•  Experimental  setup  and  testing  protocol  were  completed. 

•  Formulation  of  the  program  for  data  analysis  was  completed. 

•  Validation  of  the  program  and  model  was  completed. 

•  Determining  the  final  marker  model  for  each  specific  prosthesis  design  was  completed. 

•  Determining  a  resultant  marker  model  for  all  tested  running  prosth  esis  designs  was 
completed. 

•  This  research  determined  that  neither  the  number  of  markers  placed  nor  the 
placement  of  the  markers  on  a  running-specific  prosthesis  influence  the  force  and 
torque  transfer  estimations  to  the  proximal  end  of  the  prosthesis  during  stance 
phase. 

•  A  manuscript  for  journal  publicatio  n  is  currently  in  re  view.  A  draft  of  this  manuscript  is 
included  as  Appendix  II:  Manuscript  Draft  1  of  this  document. 
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Manuscripts  Supported  by  this  Award 

1.  Baum  BS,  Koh  K,  Linberg  A,  Tian  A,  Kim  H,  Hsieh  A,  Wolf  EJ,  Shim  JK.  (In  Review).  Marker 
placement  on  running-specific  prostheses  does  not  affect  kinetics.  Journal  of  Biomechanics. 

2.  Baum  BS,  Schultz  MP,  Tian  A,  Shefter  B,  Shim  JK.  (In  Preparation).  Determining  the  inertial 
properties  of  running-specific  prostheses.  Archives  of  Physical  Medicine  and  Rehabilitation. 


Meeting  Abstracts  Supported  by  this  Award 

1.  Baum  BS,  Shim  JK.  (2009)  Optimization  and  validation  of  a  biomechanical  model  for 
running-specific  prostheses.  Research  Interaction  Day,  University  of  Maryland,  College 
Park  (September  18,  College  Park,  MD). 

2.  Baum  BS,  Borjian  R,  Kim  YS,  Linb  erg  A,  Shim  JK.  (2010)  Optimization  and  validation  of  a 
biomechanical  model  for  analyzing  running-specific  prostheses.  The  26  th  Southern 
Biomedical  Engineering  Conference  (April  30-May  2,  College  Park,  MD). 

3.  Baum  BS,  Borjian  R,  Linberg  A,  Koh  K,  Shim  JK.  (2011)  Optimization  and  validation  of  a 
biomechanical  model  for  running-specific  prostheses.  The  15th  Annual  Meeting  of  the  Gait 
and  Clinical  Movement  Analysis  Society  (April  26-29,  Bethesda,  MD). 
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Joint  Control  Adaptations  during  Running  in  Individuals  with  Lower  Extremity  Amputation, 
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1 .  Brian  S.  Baum,  Doctor  of  Philosophy  degree,  Anticipated  March  2012.  Dissertation:  Kinetics 
in  Individuals  with  Unilateral  Transtibial  Amputations  using  Running-Specific  Prostheses. 
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Conclusions 


The  research  project  has  been  completed  as  proposed  and  a  manuscript  with  the  final  data  is  in 
review.  Some  unexpected  difficulties  (e.g.  procurement  of  prostheses;  detailed  in  the  August 
2010  Annual  Report  for  this  project  )  delayed  portions  of  the  research  in  the  first  year  of  th  e 
project.  A  no-cost  extension  was  granted  based  on  these  issues.  The  issues  were  resolved  and 
did  not  affect  the  overall  successful  completion  of  the  project. 

The  data  shown  in  the  Body  section  allowed  us  to  identify  the  iterations  (combinations  of  marker 
placements)  that  yield  ed  acceptable  error  b  etween  the  estimated  (via  inverse  dynamics 
calculations)  and  dire  ctly  measured  (via  load  cell)  force  and  moment  values.  These  data 
indicate  that  marker  placements  on  running-specific  prosth  eses  can  be  flexible  if  only  kinetic 
analyses  are  desired;  however,  more  specific  marker  placements  are  recommended  if  kinematic 
information  about  the  prosthesis  compression  is  of  interest.  These  data  will  guide  future 
research  and  provide  greater  confidence  in  reported  kinetic  results  in  the  past  literature. 

The  data  indicate  that  the  marker  combinations  tested  result  in  errors  of  less  than  1 .9%  for  force 
and  moment  calculations  for  all  running  prosth  esis  designs.  These  data  suggest  t  hat  placing 
one  marker  on  a  runnin  g-specific  prosthesis  is  sufficient  for  accurate  joint  kinetic  analyses  for 
the  proximal  joints  during  stance  phase  of  running.  Additional  resear  ch  investigating  stance 
phase  joint  kinetics  during  amputee  running  may  now  be  conducted  with  the  knowledge  that 
marker  placement  does  not  greatly  influence  these  outcome  measures. 

This  knowledge  also  allows  a  la  rger  number  of  resear  ch  laboratories  to  perf  orm  running 
analyses  since  fewer  motion  capture  camera  s  would  be  needed  to  perform  the  a  nalysis.  This 
will  also  dramatically  reduce  the  setup  time  (fewer  markers  =  less  t  ime  spent  during  setup)  and 
the  impact  on  the  individual  being  tested. 

Please  see  Appendix  II:  Manuscript  Draft  1  for  a  more  detailed  report  of  the  study  conclusions. 

The  current  study  used  geometric  estimations  to  predict  the  inertial  properties  of  the  prosthesis 
subsegments.  The  inert  ial  properties  of  differe  nt  running-specific  pro  sthesis  models  has  not 
been  reported  in  the  literature,  yet  these  data  are  necessary  as  inputs  to  the  inver  se  dynamics 
equations  of  motion.  Additional  re  search  is  w  arranted  to  accurately  determine  these  inertia  I 
parameters  of  running-specific  prostheses  to  be  used  in  conjunction  with  the  marker  placement 
data  reported  in  the  current  study  for  the  analysis  of  overground  running  kinetics  in  the  amputee 
population.  As  an  extension  of  the  current  study,  we  performed  inertial  testing  of  the  running- 
specific  prostheses  and  a  draft  of  a  manuscript  resulting  from  this  work  is  included  as  Appendix 
II:  Manuscript  Draft  2.  These  da  ta  will  provi  de  researchers  with  g  uidelines  to  assist  with 
determining  the  inertial  properties  o  f  running-specific  prostheses  or  with  estimations  to  use  in 
their  own  analyses  if  they  cannot  perform  inertial  testing  at  their  facilities. 
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Marker  Placement  on  Running-Specific  Prostheses  Does  Not  Affect  Kinetics 

Abstract 

Gait  analyses  for  individuals  wearing  running-specific  prostheses  (RSPs)  are  currently 
performed  by  placing  reflective  markers  arbitrarily  on  the  RSP  and  inverse  dynamics  techniques 
are  then  used  to  estimate  joint  kinetic  data.  Marker  placements  on  RSPs  have  not  been  validated 
for  accuracy  in  estimating  joint  kinetic  data  and  potential  errors  within  these  estimations  are 
unknown.  This  study  examined  how  varying  marker  placements  on  RSPs  affect  proximal  kinetic 
estimations  during  an  axial  loading  task.  Reflective  markers  were  placed  every  2  cm  along  the 
lateral  aspects  of  four  RSP  models  with  three  different  stiffness  categories  each  (12  total  RSPs). 
Prostheses  were  neutrally  aligned  in  a  material  testing  system  (MTS)  between  two  load  cells  that 
measured  proximal  applied  forces  and  moments  and  distal  ground  reaction  forces.  RSPs  were 
axially  loaded  to  2500N  of  force  to  simulate  peak  running  loads.  Inverse  dynamics  estimated 
force  transfers  from  the  ground  to  the  proximal  endpoint  of  the  prostheses  through  the  segments 
defined  by  reflective  markers.  Errors  between  estimated  and  applied  values  at  the  proximal 
endpoint  were  calculated  for  every  combination  of  markers.  Regardless  of  the  number  of  markers 
or  their  placement  on  the  RSPs,  joint  kinetic  estimations  resulted  in  root  mean  square  errors  less 
than  10  N  (1%),  17.5  N  (0.75%),  and  2.5  Nm  (1.6%)  for  AP  force,  vertical  force,  and  flexion 
moment,  respectively  as  compared  to  the  directly  measured  values  at  the  proximal  end  of  the 
prostheses.  The  results  suggest  that  placing  a  single  marker  on  an  RSP  is  sufficient  for  accurate 
stance  phase  kinetic  analyses. 


Introduction 


During  three-dimensional  gait  analyses,  reflective  markers  are  placed  on  anatomical 
landmarks  to  estimate  the  positions  of  joint  centers  and  to  define  the  body  segment  motions.  The 
distal  joint  motion  data  along  with  ground  reaction  force  data  from  a  force  platform  can  be  used 
as  inputs  to  inverse  dynamics  equations  to  estimate  proximal  joint  kinetic  values.  In  locomotion 
studies  using  prostheses,  markers  defining  the  most  distal  joint  axis,  usually  the  ankle,  are 
generally  affixed  to  spots  on  the  prosthetic  foot  that  mimic  the  relative  marker  location  on  the 
intact  foot  and  ankle  complex1'6.  Prostheses  are  often  modeled  anthropometrically  like  an  intact 
limb  even  though  these  devices  may  not  have  the  same  architecture  or  landmarks  '  . 

With  the  development  of  running-specific  prostheses  (RSPs),  new  prosthetic  foot  designs 
have  emerged  that  no  longer  resemble  the  human  foot.  Many  of  the  designs  resemble  a  “C”  or 
“L”  shape  at  the  distal  end  of  the  limb,  which  allows  the  prosthesis  to  flex  and  return  more 
energy  for  propulsion  during  running,  similar  to  a  spring.  Placing  multiple  markers  to  model 
RSPs  as  multisegmented  objects  during  amputee  locomotion  studies  provides  a  great  challenge 
since  definitive  joint  axes  may  not  exist  within  the  prosthetic  foot  design,  yet  modeling  RSPs  as 
single  rigid  objects  may  not  be  appropriate  since  these  devices  are  designed  to  flex  throughout 
their  length.  In  the  face  of  these  challenges,  many  researchers  analyze  these  prostheses  using 
similar  biomechanical  analysis  methods  as  have  been  employed  in  prosthetic  feet  designed  for 
walking  and  intact  feet.  Studies  investigating  running  with  RSPs  have  estimated  the  prosthetic 
limb  “ankle”  joint  to  be  either  at  the  same  relative  position  as  the  intact  limb’s  ankle  joint  (Figure 
la)  or  the  most  acute  point  on  the  prosthesis  curvature  (i.e.,  the  greatest  curvature;  Figure  lb)6, 10, 
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***  INSERT  FIGURE  1  ABOUT  HERE  *** 


These  estimations  have  not  been  validated  and  potentially  result  in  significant  errors 
within  the  kinetic  calculations  and  subsequent  interpretations  of  results.  Using  the  intact  limb  as 
a  reference  for  marker  placement  also  excludes  such  a  model  from  use  on  individuals  with 
bilateral  amputations.  Consequently,  researchers  need  to  know  how  marker  placement  on  RSPs 
affects  proximal  joint  kinetic  estimations  so  models  can  be  created  for  use  with  different  RSP 
designs  and  can  be  utilized  in  individuals  with  bilateral  lower  extremity  amputations.  An 
accurate  model  will  provide  data  that  can  be  interpreted  with  confidence  and  is  needed  to 
produce  biomechanical  and  physiological  data  necessary  to  identify  optimal  running  techniques, 
prosthetic  alignment,  prosthetic  designs,  training  regimens,  and  energy  efficiency. 

In  addition  to  the  different  designs  of  running-specific  prostheses,  each  of  these  devices 
are  manufactured  in  different  stiffness  categories  that  are  generally  prescribed  based  on  an 
individual’s  body  weight  and  general  activity  level.  A  heavier  person  is  typically  prescribed  a 
RSP  with  a  higher  category  of  stiffness  (higher  categories  correspond  to  greater  prosthesis 
stiffness).  Studies  investigating  prosthesis  stiffness  indicate  that  the  stiffness  affects 
performance,  so  body  weight  and  general  activity  level  may  be  insufficient  guidelines  for 
prescribing  a  stiffness  category.  A  stiffer  forefoot,  wider  c-curve,  and  thinner  lay-up  resulted  in 
individuals  with  lower  extremity  amputation  (ILEA)  running  their  fastest  sprint  times12,  which 
suggests  that  sprint  speed  can  be  a  function  of  stiffness  and  prosthetic  foot  shape13.  Using  a 
greater  category  of  stiffness  may  also  improve  gait  symmetry  values  for  ILEA  with  transtibial 
amputation14,  15,  but  it  has  also  been  shown  to  reduce  energy  efficiency16.  These  data  suggest 
that  different  prosthesis  stiffness  categories  could  affect  the  performance  of  the  prosthesis  and 
therefore  the  force  and  moment  transfer  through  the  device. 


The  aim  of  this  study  was  to  investigate  the  effect  of  marker  placement  on  errors  between 
calculated  force  and  moment  values  from  inverse  dynamics  estimations  and  directly  measured 
values  for  specific  running  prosthesis  designs  and  stiffness  categories.  We  hypothesized  that 
more  markers  placed  on  any  running-specific  prosthesis  with  any  stiffness  category  would  result 
in  the  least  error  in  proximal  force  and  moment  estimations. 

Methods 

Four  of  the  most  commonly  prescribed  running-specific  prostheses  currently  available  on 
the  market  were  tested  including  the  1E90  Sprinter  (OttoBock  Inc.),  Flex-Run  (Ossur),  Cheetah 
(Ossur)  and  Nitro  Running  Foot  (Freedom  Innovations)  (Figure  2).  Three  different  stiffness 
categories  were  also  tested  for  each  prosthetic  design  to  identify  whether  prosthetic  stiffness 
affects  marker  placement  results.  Stiffness  categories,  presented  in  Table  1,  were  chosen  to 
reflect  a  common  range  of  stiffnesses  that  might  be  prescribed.  OttoBock  does  not  use  the  term 
“category”  to  reflect  stiffness,  rather  different  prosthesis  stiffnesses  are  reflected  by  the  target 
weight  of  the  person  using  the  device.  Prostheses  were  aligned  neutrally  according  to 
manufacturer  recommendations  with  their  proximal  ends  attached  to  a  six-degree-of-freedom 
load  cell  (Bertec  PY6,  Columbus,  OH)  that  was  connected  to  the  arm  of  a  material  testing  system 
(MTS,  Eden  Prairie,  MN).  A  second  load  cell  was  secured  to  the  base  of  the  MTS  (Figure  3). 
The  prostheses  were  cyclically  loaded  for  ten  cycles  with  axial  forces  up  to  2,500  N  to  simulate 
peak  vertical  forces  commonly  observed  during  running17'19  (approximately  three  times  the  body 
weight  of  a  75  kg  person).  The  load  cells  sampled  data  at  1000  Hz  and  measured  the  forces  and 
moments  at  the  point  of  load  application  proximal  to  the  prostheses  (applied  load)  and  the 
reaction  forces  distal  to  the  prostheses  (ground  reaction  forces). 


***  INSERT  TABLE  1  ABOUT  HERE  *** 


***INSERT  FIGURE  2  ABOUT  HERE  *** 


Reflective  markers  were  placed  at  2  cm  intervals  along  the  lateral  aspect  of  the  keel  of 
each  running-specific  prosthesis  (see  Figure  3).  Reflective  markers  were  also  placed 
orthogonally  on  the  anterior,  lateral,  and  medial  aspect  of  the  “head”  of  the  prosthesis,  at  the 
point  of  connection  to  the  socket  or  pylon,  in  order  to  define  the  local  coordinate  system  of  the 
prosthesis.  Three  additional  markers  were  placed  along  the  midline  of  each  prosthesis  to  define  a 
plane  to  which  the  keel  markers  were  projected  for  further  analysis.  An  8-camera  motion  capture 
system  (Vicon,  Oxford,  UK)  with  a  capture  frequency  of  500  Hz  was  used  to  collect  the  3-D 
positional  data  of  the  markers  during  each  trial.  Two  consecutive  projected  midline  markers 
defined  individual  segments  of  the  prosthesis  (assumed  to  be  rigid)  and  consecutive  segments 
shared  a  common  marker.  The  joint  between  these  segments  was  assumed  as  a  hinge  joint. 
Standard  inverse  dynamics  calculations  were  made  to  estimate  the  force  and  torque  transfer 
from  the  ground  reaction  force,  through  the  defined  prosthesis  segments,  and  to  the  point  of  load 
application  proximal  to  the  prosthesis. 

***  INSERT  FIGURE  3  ABOUT  HERE  *** 

***  INSERT  FIGURE  4  ABOUT  HERE  *** 

Prosthesis  thickness  was  measured  at  each  marker  position  using  digital  calipers,  and 
prosthesis  width  at  each  position  was  calculated  as  twice  the  distance  between  the  marker  and  its 
midline  projection.  Prosthesis  segments  were  defined  by  two  consecutive  markers  and  were 
considered  as  rigid  trapezoidal  cuboids  (see  Figure  4).  The  center  of  mass  along  the  width  and 
thickness  of  each  segment  were  determined  from  half  the  average  width  and  thickness, 
respectively.  The  center  of  mass  position  along  the  long  axis  (length)  of  each  segment  was 
determined  by  equation  1 : 
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where  Wd  and  wp  are  the  distal  and  proximal  end  widths  and  /  is  the  segment  length. 

Segment  volumes  were  estimated  as  a  trapezoidal  cuboid  volume  and  the  total  volume  of 
each  prosthesis  was  estimated  by  summing  all  segment  volumes.  Mass  was  assumed  to  be  evenly 
distributed  throughout  each  RSP  such  that  the  ratio  of  segment  volume  to  total  volume  equaled 
the  ratio  of  segment  mass  to  total  mass,  and  segment  masses  were  determined  accordingly.  The 
inertial  properties  of  each  prosthesis  segment  were  estimated  using  assumptions  based  on  a 
trapezoidal  cuboid.  Each  segment  length  was  integrated  across  200  subsegments.  The  principal 
axis  moments  of  inertia  of  each  segment  were  estimated  by  equations  2-4: 
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where  nij  is  the  mass,  /,  is  the  length,  vt-v  is  the  width,  U  is  the  thickness,  r,  is  the  distance  between 
the  subsegment  center  of  mass  and  the  segment  center  of  mass  for  each  subsegment  i, 
respectively. 

The  angles  between  each  set  of  three  consecutive  markers  were  calculated  throughout  the 
cyclic  loading  and  the  range  of  angle  change  was  determined  at  each  marker  “joint”.  Markers 
representing  joints  that  had  an  angular  change  of  less  than  one  degree  were  removed  from  further 
analyses  as  they  were  considered  as  part  of  a  larger  rigid  segment.  The  remaining  markers  were 
used  for  the  model  analysis.  The  one  degree  threshold  was  determined  from  the  marker  position 
error  of  the  motion  capture  system. 


The  difference  between  force  and  moment  values  at  the  point  of  load  application  from  the 
estimated  inverse  dynamics  calculations  and  the  directly  measured  values  from  the  top  load  cell 
was  considered  model  error.  Force  and  moment  estimations  were  made  with  every  combination 
of  remaining  markers  giving  a  resultant  error  value  for  each  combination.  Error  was  calculated 
for  each  loading  cycle  using  equations  5  and  6  for  root  mean  squared  error  (RMSE)  and 
normalized  RMSE  (NRMSE),  respectively. 


RMSE  - 

ft 

[5] 

NRMSE  =  ,  - -  ft 

[6] 

where  Km  represents  the  directly  measured  kinetic  values  (force  or  moment)  from  the  upper  load 
cell,  Kc  represents  the  calculated  kinetic  values  from  inverse  dynamics  equations,  n  is  the  number 
of  data  points  in  the  loading  cycle,  and  max  and  min  represent  the  maximum  and  minimum 
values  within  the  loading  cycle,  respectively. 

Uncertainty  Analysis 

The  effects  of  input  variable  uncertainties  were  used  to  estimate  the  uncertainty  of  the 

resultant  joint  force  and  moment  variables  via  an  error  analysis  method.  The  upper  bound  of 
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uncertainty  in  the  result  ur  was  calculated  according  to  equation  7:  ’ 

= 

where  R  is  the  resultant  value  (e.g.  joint  force  or  joint  moment),  x;  is  the  ith  input  variable  in 
predicting  R,  and  Ax;  is  the  error  associated  with  input  variable  Xj.  Primary  sources  of  error  (Ax;) 
included  errors  related  to  load  cell  measurements  and  marker  noise  which  affects  segmental 
kinematic  parameters.  The  uncertainties  estimated  for  AP  force,  vertical  force,  and  flexion 
moment  for  each  prosthesis  are  shown  in  Table  2. 


***  INSERT  TABLE  2  ABOUT  HERE  *** 


Results 

Calculated  values  and  error  data  are  presented  for  anteroposterior  (AP)  forces,  vertical 
forces,  and  flexion  moments  during  the  cyclical  loading  trials  for  each  prosthesis.  Mediolateral 
(ML)  forces,  ML  rotational  moments,  and  intemal/extemal  rotational  moments  are  not  presented 
since  the  axial  loading  of  the  prostheses  produced  minimal  forces  and  moments  along  and  about 
these  axes,  respectively. 

Regardless  of  the  number  of  markers  or  their  placement  on  the  various  RSPs,  force  and 
moment  calculations  using  inverse  dynamics  techniques  resulted  in  errors  of  less  than  1.6%  as 
compared  to  the  directly  measured  values  (Table  3).  Directly  measured  and  calculated  AP  force, 
vertical  force,  and  flexion  moment  values  are  presented  in  Figure  5.  Raw  errors  between  the 
directly  measured  and  calculated  forces  and  moments  are  presented  in  Figure  6. 

***  INSERT  TABLE  3  ABOUT  HERE  *** 

***  INSERT  FIGURE  5  ABOUT  HERE  *** 

***  INSERT  FIGURE  6  ABOUT  HERE  *** 

The  Freedom  Innovations  Nitro  prosthesis  had  a  maximal  RMSE  range  of  0.26  N  (AP 
force),  4.45  N  (vertical  force),  and  1 .02  Nm  (flexion  moment)  and  a  maximal  NRMSE  range  of 
0.02%,  0.17%,  and  0.86%  for  AP  force,  vertical  force,  and  flexion  moment,  respectively  across 
all  stiffness  categories  and  all  tested  combinations  of  markers.  The  Ossur  Flex-Run  prosthesis 
had  a  maximal  RMSE  range  of  4.37  N  (AP  force),  5.88  N  (vertical  force),  and  1.05  Nm  (flexion 
moment)  and  a  maximal  NRMSE  range  of  0.37%,  0.28%,  and  0.56%  for  AP  force,  vertical  force, 
and  flexion  moment,  respectively  across  all  stiffness  categories  and  all  tested  combinations  of 


markers.  The  Ossur  Cheetah  prosthesis  had  a  maximal  RMSE  range  of  0.99  N  (AP  force),  9.38  N 
(vertical  force),  and  0.73  Nm  (flexion  moment)  and  a  maximal  NRMSE  range  of  0.12%,  0.44%, 
and  0.53%  for  AP  force,  vertical  force,  and  flexion  moment,  respectively  across  all  stiffness 
categories  and  all  tested  combinations  of  markers.  The  Ottobock  1E90  prosthesis  had  a  maximal 
RMSE  range  of  0.48  N  (AP  force),  7.54  N  (vertical  force),  and  0.54  Nm  (flexion  moment)  and  a 
maximal  NRMSE  range  of  0.07%,  0.35%,  and  0.31%  for  AP  force,  vertical  force,  and  flexion 
moment,  respectively  across  all  stiffness  categories  and  all  tested  combinations  of  markers. 

The  average  AP  force,  vertical  force,  and  flexion  moment  RMSE  values  for  each 
combination  according  to  the  number  of  markers  on  the  prostheses  are  shown  in  Figure  7.  These 
data  show  little  difference  in  RMSE  values  in  kinetic  variables  regardless  of  the  number  of 
markers  on  a  prosthesis. 

***  INSERT  FIGURE  7  ABOUT  HERE  *** 

Discussion 

This  study  examined  the  effects  of  marker  placement  on  proximal  kinetic  estimations 
using  inverse  dynamics  during  a  cyclic  loading  task  with  running-specific  prostheses.  The  results 
of  this  study  indicated  that  root  mean  square  errors  between  the  directly  measured  and  calculated 
kinetic  variables  were  less  than  18N  for  vertical  forces,  10N  for  AP  forces,  and  2Nm  for  flexion 
moments.  Considering  peak  values  of  approximately  2500N,  700N,  and  120Nm  for  vertical 
force,  AP  force,  and  flexion  moment,  respectively,  NRMSE  values  were  less  than  or  equal  to 
1.6%  for  all  combinations  of  marker  placements  across  all  prostheses  investigated.  These  low 
errors  indicate  that  using  any  combination  of  markers  would  result  in  proximal  joint  kinetic 
estimations  with  reasonable  errors  for  a  gait  analysis. 


To  investigate  whether  placing  a  particular  number  of  markers  on  a  prosthesis  would 
have  an  effect  on  the  outcome  variables,  the  average  RMSE  values  were  calculated  for  all 
combinations  of  a  particular  number  of  markers.  For  example,  errors  were  averaged  for  all 
combinations  of  one  marker  placed  on  a  prosthesis,  all  combinations  of  two  markers  placed  on  a 
prosthesis,  all  combinations  of  three  markers  placed  on  a  prosthesis,  etc.  This  data  showed 
similar  average  RMSE  values  for  the  combinations  of  particular  numbers  of  markers  placed  on 
the  prostheses  (see  Figure  7).  This  information  combined  with  the  relatively  small  range  of  errors 
across  all  tested  marker  combinations  suggest  that  the  number  and  placement  of  markers  on  any 
of  the  tested  running-specific  prostheses  does  not  greatly  influence  the  estimation  of  force  and 
moment  transfer  through  the  prostheses.  The  estimated  uncertainty  values  for  AP  force,  vertical 
force,  and  flexion  moment  were  less  than  1%  of  the  peak  force  and  moment  values.  The  range  of 
RMSE  across  marker  combinations  for  any  prosthesis  was  similar  to  the  uncertainty  values  for 
that  prosthesis  indicating  that  all  marker  combinations  had  similar  RMSE  values. 

One  possible  explanation  for  the  small  change  in  error  across  different  marker  positions 
is  that  the  magnitudes  of  the  ground  reaction  forces  during  running  are  very  large  in  comparison 
to  the  accelerations  and  inertial  properties  of  the  running-specific  prostheses.  During  force 
transfer  from  the  ground  through  the  prostheses,  the  centers  of  mass  of  the  prosthesis 
subsegments  do  not  change  dramatically  (low  accelerations).  Therefore,  the  ground  reaction 
force  transfers  nearly  unattenuated  through  the  prosthesis.  The  ground  reaction  forces  also 
generate  torques  that  account  for  nearly  all  of  the  estimated  proximal  joint  moments  while  the 
moments  of  inertia  and  the  angular  velocities  of  the  prosthetic  segments  contribute  relatively 


little  to  stance  phase  kinetics. 


These  data  suggest  that  kinetic  data  calculated  from  prior  research  with  RSPs  may  be 
interpreted  with  greater  confidence.  Placing  markers  at  the  same  relative  position  as  the  intact 
limb’s  ankle  joint  or  the  most  acute  point  on  the  prosthesis  curvature6, 10, 11  should  yield  similar 
results  in  resultant  kinetic  values  proximal  to  the  prosthesis.  However,  for  consistency  and 
flexibility  in  modeling,  it  is  recommended  that  markers  are  placed  according  to  the  prosthesis 
architecture  rather  than  intact  limb  architecture.  This  will  allow  markers  to  be  placed  on  the  same 
location  of  a  particular  prosthesis  from  subject  to  subject  and  will  allow  for  the  study  of  ILEA 
with  bilateral  amputations  and  comparison  of  these  individuals  with  those  with  unilateral 
amputations. 

Most  motion  capture  laboratories  have  a  limited  number  of  cameras  and  may  have 
difficulty  tracking  a  large  number  of  markers  placed  closely  together  during  activities  such  as 
running.  This  limits  the  number  of  markers  that  researchers  can  feasibly  place  on  the  keel  of  a 
running-specific  prosthesis,  especially  considering  that  the  thin  profile  of  such  prostheses  often 
necessitates  using  markers  with  small  diameters.  Furthermore,  motion  capture  of  overground 
running  requires  a  large  capture  volume,  and  optimizing  camera  placement  for  large  volumes 
reduces  the  effectiveness  of  these  systems  to  capture  small  markers  in  close  proximity  to  each 
other.  Utilizing  a  minimal  marker  set  for  running- specific  prostheses  will  enable  widespread  use 
of  such  a  model  regardless  of  the  number  of  cameras  available  to  a  laboratory  and  to  allow  for 
both  overground  and  treadmill  data  collections  while  using  the  same  model.  Additionally,  fewer 
markers  on  a  prosthesis  makes  setup  less  tedious  and  saves  testing  time. 

Several  limitations  exist  in  this  study.  First,  only  axial  loading  was  performed  on  the 
prostheses,  whereas  when  running,  the  prostheses  are  loaded  while  rolling  forward,  which  would 
produce  different  loading  patterns  and  potentially  different  prosthetic  bending.  This  could  affect 


the  recommended  marker  placements  on  the  prostheses.  However,  the  overall  ground  reaction 
forces  during  running  are  still  much  larger  than  the  inertial  properties  of  the  prostheses,  so  it  is 
anticipated  that  for  kinetic  analyses,  the  results  presented  in  this  study  would  generalize  to 
overground  running.  However,  due  to  the  axial  loading,  this  study  only  presented  AP  force, 
vertical  force,  and  flexion  moment  results.  Validation  of  the  marker  models  is  still  needed  for 
mediolateral  forces,  varus/valgus  moments,  and  intemal/extemal  rotational  moments.  Additional 
studies  are  warranted  to  investigate  these  kinetic  parameters  using  either  a  6-degree-of- freedom 
material  testing  system  that  could  mimic  the  prosthetic  roll-over  during  running  or  direct  load 
measurements  at  the  proximal  end  of  the  prosthesis  during  running.  An  additional  limitation  of 
this  study  is  that  only  stance  phase  loading  was  investigated.  The  inertial  effects  of  the  running 
prostheses  during  swing  phase  are  most  likely  not  trivial,  so  accurate  measures  of  mass,  center  of 
mass  position,  and  moments  of  inertia  are  needed  to  accurately  estimate  the  joint  kinetic  values 
proximal  to  the  prostheses.  Future  studies  are  needed  to  accurately  measure  and  predict  the 
inertial  properties  and  effects  of  running-specific  prostheses  during  the  running  swing  phase. 

The  development  and  validation  of  an  accurate  biomechanical  model  for  use  with 
running-specific  prostheses  allows  researchers  to  fully  examine  the  kinematic  and  kinetic 
adaptations  that  occur  during  running  in  ILEA.  Extremely  limited  information  is  available  in  the 
literature  to  guide  clinicians  in  aligning,  prescribing,  or  rehabilitating  ILEA  who  wish  to  run.  For 
example,  it  is  currently  unknown  whether  running  with  running-specific  prostheses  poses  an 
increased  risk  for  injury  in  the  residual  limb  joints  or  joints  in  the  contralateral  limb.  ILEA  are 
already  at  greater  risk  of  degenerative  joint  diseases  such  as  osteoarthritis  (OA)23,  and  the  larger 
forces  generated  during  running  could  promote  the  development  and  progression  of  these 
diseases.  Prior  research  supports  that  OA  may  initiate  in  joints  that  experience  a  traumatic  or 


24 

chronic  event  (such  as  amputation  due  to  injury  or  disease)  that  causes  kinematic  changes  .  The 
rate  of  OA  progression  is  currently  thought  to  be  associated  with  increased  loads  during 
ambulation  ’  .  Identifying  running  techniques,  prosthetic  alignments,  or  new  prosthetic  designs 
that  reduce  peak  lower  extremity  joint  loading  may  reduce  the  risk  of  developing  and  progressing 
OA. 

Additional  research  needs  include  investigating  the  effects  of  various  prosthetic 
components  in  meeting  different  running  goals,  and  determining  optimal  prosthetic  alignment  so 
as  to  minimize  asymmetries  and  maximize  energy  efficiency  during  running. 

Conclusions 

Regardless  of  the  number  of  markers  or  their  placement  on  the  various  RSPs,  joint  kinetic 
estimations  resulted  in  root  mean  square  errors  less  than  10  N  (1%),  17.5  N  (0.75%),  and  2.5  Nm 
(1.6%)  for  AP  force,  vertical  force,  and  flexion  moment,  respectively  as  compared  to  the  directly 
measured  values  at  the  proximal  end  of  the  prostheses  during  axial  loading  up  to  2,500  N.  This 
affords  researchers  the  flexibility  to  place  markers  conveniently  on  running-specific  prostheses 
and  still  confidently  estimate  joint  kinetic  data  during  the  stance  phase  of  running. 
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Table  1.  Stiffness  categories  used  for  each  prosthesis  during  testing  and  the 
manufacturer  recommended  body  mass  range  associated  with  each  category. 

Prosthesis  Stiffness  Category  (body  mass  range) 


Flex-Run 

3  (53-59  kg) 

5  (69-77  kg) 

7  (89-100  kg) 

Nitro 

3  (60-68  kg) 

6  (89-100  kg) 

7  (101-1 16  kg) 

Cheetah 

3  (60-68  kg) 

5  (78-88  kg) 

7  (101-1 16  kg) 

1E90 

140  lb  (63.6  kg) 

185  lb  (84.1kg) 

235  lb  (106.8  kg) 

Table  3.  Error  ranges  (minimum  to  maximum  RMSE  and  NRMSE)  of  all  combinations  of  markers  for  the  estimated  kinetic  values  from 


inverse  dynamics  equations. 

Freedom  Innovations 

Ossur 

Ossur 

Ottobock 

Nitro 

Flex-Run 

Cheetah 

1E90 

Stiffness  Category: 

Cat  3 

Cat  6 

Cat  7 

Cat  3 

Cat  5 

Cat  7 

Cat  3 

Cat  5 

Cat  7 

1401b 

185  lb 

235  1b 

RMSE 

6.78- 

4.32- 

5.66- 

5.17- 

9.39- 

3.23- 

2.17- 

1.36- 

2.46- 

0.29- 

6.27- 

5.28- 

AP 

(N) 

6.92 

4.50 

5.92 

9.54 

9.54 

4.21 

2.27 

1.59 

3.45 

0.77 

6.46 

5.56 

Force 

NRMSE 

0.68- 

0.61- 

0.54- 

0.56- 

0.92- 

0.35- 

0.31- 

0.34- 

0.29- 

0.05- 

1.43- 

0.90- 

(%) 

0.69 

0.63 

0.56 

0.93 

0.93 

0.45 

0.32 

0.40 

0.41 

0.12 

1.47 

0.95 

RMSE 

16.37- 

11.55- 

14.39- 

11.05- 

10.19- 

11.85- 

10.27- 

7.57- 

9.22- 

6.88- 

7.17- 

7.16- 

Vertical 

(N) 

16.80 

16.00 

15.85 

16.93 

11.95 

14.18 

17.41 

16.95 

17.49 

11.93 

9.61 

14.70 

Force 

NRMSE 

0.64- 

0.45- 

0.50- 

0.44- 

0.41- 

0.49- 

0.30- 

0.36- 

0.28- 

0.41- 

0.47- 

0.33- 

(%) 

0.66 

0.62 

0.55 

0.72 

0.48 

0.59 

0.50 

0.80 

0.53 

0.71 

0.63 

0.68 

RMSE 

0.98- 

0.81- 

1.27- 

0.63- 

0.78- 

0.89- 

0.91- 

0.87- 

0.66- 

2.02- 

0.98- 

1.03- 

Flexion 

(Nm) 

1.36 

1.83 

1.99 

1.59 

1.09 

1.94 

1.13 

1.14 

1.39 

2.32 

1.52 

1.38 

Moment 

NRMSE 

0.52- 

0.67- 

0.72- 

0.31- 

0.37- 

0.48- 

0.67- 

0.53- 

0.46- 

0.75- 

0.58- 

0.38- 

(%) 

0.71 

1.53 

1.14 

0.78 

0.51 

1.04 

0.83 

0.70 

0.99 

0.86 

0.89 

0.51 

*Notes:  RMSE  =  root  mean  square  error,  NRMSE  =  normalized  root  mean  square  error 
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Figure  1.  Literature  has  reported  marker  placement  for  running  prostheses  at  (a)  the  height  of 
the  intact  limb’s  lateral  malleolus  or  (b)  the  point  at  which  the  radius  of  the  prosthesis  is  most 
acute. 


a.  Freedom  Innovations  b.  Ossur  Flex-Run  c.  Ossur  Cheetah  d.  Ottobock  Sprinter 

Nitro  1E90 

Figure  2.  Prostheses  used  for  mechanical  testing. 


Figure  3.  Marker  placement  on  a  running-specific  prosthesis  and  its  position  in  an  MTS 
machine  between  two  load  cells.  Fewer  markers  than  actual  are  shown  in  the  illustration  for 
clarity.  The  red  dot  indicates  the  point  of  load  application,  measured  by  the  upper  load  cell. 
The  lower  load  cell  measured  ground  reaction  force  (GRF).  The  red  arrows  represent  the  input 
and  GRF  force  vectors. 


td 


Figure  4.  Schematic  of  segment  definitions  within  each  prosthesis.  Blue  circles  represent 
markers,  the  red  circle  represents  the  segment  center  of  mass.  The  axis  defines  the  segment 
local  coordinate  system  with  its  origin  at  the  center  of  mass.  Segment  length  (/)  is  also  shown 
along  with  the  width  (w)  and  thickness  (t)  at  the  proximal  (p)  and  distal  (d)  ends. 


AP  Force  Vertical  Force 


Flexion  Moment 


Figure  5.  Anteroposterior  (AP)  force,  vertical  force,  and  flexion  moment  curves  for  cyclical 
loading.  Thick  blue  lines  represent  the  directly  measured  values  from  the  upper  load  cell.  Thin 
lines  overlaid  on  the  curves  (showing  nearly  identical  patterns)  represent  calculated  values 
from  each  different  combination  of  markers.  Exemplar  data  are  from  the  Flex-Run  category  3 
prosthesis.  Other  tested  prostheses  and  stiffness  categories  showed  similar  results. 


Figure  6.  Average  anteroposterior  (AP)  force,  vertical  force,  and  flexion  moment  error  curves 
for  the  loading  cycle.  Each  curve  represents  the  difference  between  the  directly  measured 
values  from  the  upper  load  cell  and  calculated  values  from  each  combination  of  markers  on  the 
prosthesis.  Exemplar  data  are  from  the  Flex-Run  category  3  prosthesis.  Other  tested  prostheses 
and  stiffness  categories  showed  similar  results. 
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Figure  7.  Average  anteroposterior  (AP)  force,  vertical  force,  and  flexion  moment  root  mean  square  error  (RMSE)  for  each  prosthesis  across 
the  number  of  markers  on  the  prosthesis.  All  tested  combinations  with  the  number  of  markers  indicated  were  averaged  to  generate  each  data 
point.  Error  bars  represent  ±1  standard  deviation  of  all  marker  combinations  tested  for  the  number  of  markers  shown. 


Appendix  III:  Manuscript  Draft  2 


Determining  the  Inertial  Properties  of  Running-Specific  Prostheses 

Abstract 

Body  segment  inertial  properties  (mass,  center  of  mass  position,  and  moment  of  inertia) 
affect  the  joint  kinetics  and  subsequent  limb  control  during  ambulation.  Accurate  segment 
inertial  property  estimations  are  necessary  for  properly  calculating  kinetic  parameters  and 
appropriately  interpreting  these  data.  Limited  information  currently  exists  related  to  inertial 
properties  of  various  running-specific  prosthetic  designs  which  limits  the  joint  kinetic 
information  that  can  be  generated  and  interpreted  during  running  analyses.  The  purpose  of  this 
study  was  to  determine  the  inertial  properties  for  four  of  the  most  common  running-specific 
prostheses  currently  available  on  the  market  across  three  stiffness  categories  for  each  model.  A 
predictive  equation  was  developed  to  estimate  the  center  of  mass  (COM)  positions  of  the  various 
prosthetic  designs,  and  a  trifilar  pendulum  is  described  for  measuring  the  prostheses  moments  of 
inertia.  The  predictive  equation  matched  the  measured  COM  exactly  when  prosthesis-specific 
inputs  were  used  and  produced  errors  between  0. 1-2.8  cm  when  using  average  input  values 
across  prostheses  for  a  specific  design.  The  trifilar  pendulum  estimated  inertial  properties  within 
-6.21  xl0~5  kg-m2  (<  1%)  of  a  block  with  known  inertia.  Misalignments  of  the  block’s  COM  with 
the  pendulum’s  center  of  rotation  <  5cm  yielded  errors  <  0.0027  kg-m2.  Inertial  properties  about 
any  axis  varied  <  0.0038  kg-m2  within  the  tested  prosthetic  designs  due  to  different  stiffness 
categories,  although  inertial  properties  differed  more  substantially  between  different  designs. 
Inertial  estimation  errors  from  pendulum  measurements  were  less  than  or  equal  to  errors 
associated  with  various  methods  for  predicting  intact  limb  inertial  properties  indicating  that  the 


methods  and  values  presented  are  within  currently  accepted  tolerances  for  inertial  property 
estimations  for  gait  studies. 


Introduction 

Joint  kinetic  data  (e.g.  forces,  moments,  powers,  etc.)  provide  insights  into  how 
individuals  ambulate  and  how  injured  or  diseased  individuals  adapt  their  joint  control  to  achieve 
the  same  movement  goals.  Link-segment  models  allow  the  estimation  of  proximal  joint  kinetics 
from  either  distal  link  joint  kinetics  or  from  measured  distal  force  data  such  as  ground  reaction 
forces.  This  is  achieved  through  knowledge  of  the  link  motion,  inertial  properties  of  link 
segments,  and  distal  joint  kinetic  data.  Accurate  kinetic  estimations  using  link-segment  models 
depend  on  accurate  segment  inertial  property  estimations  including  mass,  center  of  mass 
position,  and  moments  of  inertia1. 

Intact  limb  inertial  properties  and  regression  equations  have  been  established  through 
cadaveric  studies2"4  and  body  scanning  methods5"7,  and  are  commonly  used  during  gait  analyses. 
However,  in  individuals  with  lower  extremity  amputation  (ILEA),  prosthetic  components  replace 
the  intact  limbs,  and  the  inertial  properties  of  the  resultant  limb-prosthesis  combination  are 
altered.  It  is  common  practice  for  researchers  to  approximate  the  inertial  properties  of  walking 
prostheses  by  using  the  values  of  intact  limbs,  since  the  prostheses  assume  similar  shapes  to  the 
foot  and  have  masses  comparable  to  those  of  the  intact  limbs.  Some  researchers  suggest  that 
modeling  prosthetic  feet  using  the  same  marker  placements  and  inertial  properties  as  intact  limbs 
produces  reasonably  acceptable  error  levels  in  gait  parameters  during  stance  phase8"10.  However, 
other  research  supports  the  notion  that  the  inertial  properties  of  prosthetic  feet  significantly 
impact  the  resultant  joint  kinematic  and  kinetic  estimations11"14  and  that  more  accurate 
estimations  of  prosthetic  inertial  properties  are  warranted.  With  the  advent  of  running-specific 


prostheses  (RSPs),  these  prosthetic  components  no  longer  resemble  the  intact  foot  and  ankle 
complex,  and  they  have  much  smaller  masses  than  the  body  parts  they  replace.  It  is  therefore 
reasonable  to  assume  that  RSPs  also  have  substantially  different  moments  of  inertia  than  intact 
limbs.  Currently  there  is  very  limited  information  on  the  inertial  properties  of  RSPs, 
measurements  of  these  properties,  and  the  effects  that  these  properties  have  on  joint 
biomechanics. 

To  our  knowledge,  only  one  study  to  date15  has  reported  any  inertial  property  values  for 
running  specific  prostheses.  This  lack  of  information  makes  it  difficult  to  compare  prosthetic 
properties  between  studies  and  to  gain  a  broader  view  of  population-wide  prosthetic  inertial 
parameters.  Multiple  methods  of  measuring  the  inertial  properties  of  objects  exist.  Genta  and 
Delprete  examined  these  methods  and  broadly  categorized  them  into  oscillitory  and 
acceleratory16.  They  reported  that  acceleratory  methods  are  more  affected  by  the  presence  of 
damping  than  are  oscillitory  methods,  and  that  torsional  and  multifilar  pendulums  are  the  most 
accurate,  capable  of  errors  less  than  1%.  Physical  pendulums  are  most  commonly  used  in  the 
prosthetic  literature  to  measure  prosthesis  moments  of  inertia.  These  pendulums  generally  rely 
on  a  joint  or  bearing  that  is  assumed  to  be  frictionless  to  make  accurate  moment  of  inertia 
measurements17, 18.  In  practice,  however,  friction  in  this  bearing  does  exist  and,  along  with  air 
resistance,  will  slow  the  period  of  oscillation  and  impart  error  in  the  inertial  estimations.  One 
solution  to  this  problem  is  to  perform  multiple  oscillation  trials  with  the  pendulum  and  use  only 
the  first  period  under  the  assumption  that  it  best  represents  the  true  period  of  oscillation17, 19.  An 
alternative  solution  is  to  use  other  pendulum  designs  that  do  not  rely  on  bearings,  such  as 
multifilar  pendulums.  A  trifilar  pendulum  is  a  form  of  multifilar  pendulum  that  utilizes  a  frame 
or  platform  suspended  from  three  equidistant  wires  about  which  rotation  occurs.  The  wire 


suspension  virtually  eliminates  the  issues  caused  by  friction  from  a  bearing  and  allows  for 
accurate  measures  of  an  object’s  moment  of  inertia  about  a  particular  axis  of  rotation.  It  also 
allows  for  measurements  over  many  more  periods  of  oscillation  than  physical  pendulums. 

The  purpose  of  this  study  was  to  determine  the  inertial  properties  (mass,  center  of  mass, 
and  principal  axis  moments  of  inertia)  for  four  of  the  most  common  running-specific  prostheses 
(RSPs)  currently  available  on  the  market.  The  construction  and  use  of  a  trifilar  pendulum  to 
estimate  these  values  are  described  in  detail.  A  secondary  aim  of  this  project  is  to  determine  the 
error  in  inertia  estimations  due  to  linear  and  rotational  malalignment  between  an  object’s  (i.e.  a 
RSP)  and  pendulum’s  center  of  mass  and  axis  of  rotation,  respectively. 

Methods 

Mass,  center  of  mass  positions,  and  principal  axis  moments  of  inertia  were  estimated 
from  four  commonly  available  running-specific  prostheses.  The  tested  models  included  the 
Freedom  Innovations  Nitro,  Ossur  Cheetah,  Ossur  Flex-Run,  and  Otto  Bock  1E90.  Three 
different  stiffness  categories  were  investigated  from  each  prosthesis  model  to  identify  whether 
inertial  differences  exist  within  each  model  type.  Categories  are  named  arbitrarily  by  the 
manufacturers,  with  increasing  stiffness  correlating  with  increasing  category  number.  The 
stiffness  categories  are  typically  prescribed  according  to  patient  body  mass  and  activity  level, 
with  greater  categories  corresponding  to  greater  weight  and  activity  intensity  (i.e.  a  greater 
stiffness  category  would  be  prescribed  for  sprinting  than  for  jogging).  Stiffness  categories  are 
also  not  standardized  across  manufacturers,  so  a  comparable  range  of  stiffness  categories  for 
each  prosthesis  was  tested.  Tested  stiffness  categories  and  their  corresponding  recommended 
body  mass  ranges  are  presented  in  Table  1. 


*****  INSERT  TABLE  1  ABOUT  HERE  ***** 


Mass  and  Center  of  Mass 

Prosthesis  masses  were  measured  using  a  standard  laboratory  scale  with  a  resolution  of 
lg.  The  center  of  mass  for  each  prosthesis  in  the  sagittal  plane  (x-z  plane)  was  measured  using  a 
reaction  board  method.  The  mass  was  assumed  to  be  evenly  distributed  throughout  the  frontal 
plane  of  the  prosthesis  such  that  the  local  y-coordinate  position  of  the  center  of  mass  would  be 
zero  for  each  prosthesis.  Based  on  these  data,  an  equation  was  developed  to  estimate  the  center 
of  mass  position  for  each  RSP  model  relative  to  the  most  proximal  (“head”)  and  most  distal 
(“toe”)  point  on  the  prosthesis.  This  equation  can  be  used  to  estimate  the  center  of  mass  of  a 
prosthesis  in  the  absence  of  a  reaction  board  or  other  equipment  needed  to  directly  measure  a 
prosthesis’  center  of  mass.  The  center  of  mass  (COM)  relates  to  the  head  and  toe  via  Equation  1 : 
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where  r  is  the  ratio  of  the  head-COM  and  head-toe  vector  magnitudes  calculated  for  the  different 
categories,  |T|  is  the  magnitude  of  the  head-toe  vector,  0  is  the  angle  between  the  head-toe  and 

head-COM  vectors  and  is  used  to  define  the  rotation  matrix,  and  up  describes  the  components  (x, 
y,  and  z)  of  the  head-toe  unit  vector.  Figure  1  shows  a  schematic  of  the  relationship  between  the 
center  of  mass,  head,  and  toe  of  a  prosthesis  with  respect  to  equation  1 . 

*****  INSERT  FIGURE  1  ABOUT  HERE  ***** 

Moments  of  Inertia 

The  inertial  properties  of  all  prostheses  were  estimated  by  placing  them  on  a  trifilar 
pendulum  (see  Figure  2)  and  measuring  the  periods  of  oscillation,  as  described  by  du  Bois  et  al  . 


The  pendulum  consisted  of  an  equilateral  plexiglass  triangle  suspended  from  its  corners  by 
equidistant  wires.  A  custom-built  aluminum  frame  served  to  support  the  trifilar  pendulum. 
Attention  was  given  to  construct  equidistant  measurements  of  the  pendulum  frame  and  to  ensure 
a  level  triangular  platform.  Frame  and  wire  contact  points  were  constructed  such  as  to  minimize 
friction  points.  A  laser  sensor  (model  BJN50-NDT,  Autonics,  Mundelein,  IL;  response  time  < 
lms)  measured  the  pendulum’s  period  of  oscillation.  The  sensor  was  aligned  such  that  the  laser 
would  be  interrupted  when  a  comer  of  the  trifilar  plate  passed  the  sensor.  Two  consecutive 
passes  determined  one  full  period  of  oscillation  of  the  pendulum. 

*****  INSERT  FIGURE  2  ABOUT  HERE  ***** 

The  triangular  platform’s  COM  was  determined  as  the  centroid  of  the  triangle 
(intersection  of  the  lines  formed  by  the  median  of  each  length  to  the  opposing  comer).  The 
center  of  mass  in  the  sagittal  plane  (the  local  x-z  plane)  for  each  prosthesis  was  estimated  by 
appropriately  orienting  and  placing  the  device  on  a  force  plate.  The  resulting  center  of  pressure 
provided  the  coordinate  of  the  COM  projection  onto  the  force  plate.  The  local  x-  and  z- 
coordinates  relative  to  the  most  proximal  end  of  the  prosthesis  (“head”)  were  then  calculated 
from  this  projection.  The  carbon  fiber  of  the  prostheses  was  assumed  to  have  evenly  distributed 
mass,  so  the  y-coordinate  of  each  prosthesis  COM  was  calculated  to  be  aligned  with  the  midline 
of  the  prosthesis. 

The  local  coordinate  system  of  each  prosthesis  was  defined  to  originate  at  the  COM 
(Figure  3).  Each  prosthesis  design  contained  a  linear  “arm”  section  at  its  proximal  end  that  was 
used  as  a  reference  to  define  the  primary  axes.  For  the  Nitro  and  Flex-Run  models,  the  most 
proximal  linear  segment  was  horizontal  when  held  in  or  near  a  neutral  alignment,  and  was  thus 
used  to  define  the  x-axis  (antero-posterior).  For  the  Cheetah  and  1E90  models,  the  proximal  end 


was  vertical,  and  was  thus  used  to  define  the  z-axis  (superior-inferior).  In  all  prostheses,  the  y- 
axis  was  parallel  to  the  width  of  the  prosthesis,  and  orthogonal  to  the  previously  defined  axis. 
The  final  principal  axis  was  defined  as  orthogonal  to  both  existing  axes. 

*****  INSERT  FIGURE  3  ABOUT  HERE  ***** 

*****  INSERT  FIGURE  4  ABOUT  HERE  ***** 


The  principal  axis  moments  of  inertia  for  each  prosthesis  were  estimated  by  measuring 
the  period  of  oscillation  from  the  trifilar  pendulum.  Prostheses  were  placed  on  the  triangular 
platform,  and  the  two  centers  of  mass  were  aligned.  In  addition,  the  prosthesis’  principal  axis 
was  aligned  with  vertical  (Figure  4).  The  platform-prosthesis  system  was  then  oscillated  about 
each  primary  axis  of  rotation.  The  moment  of  inertia  about  each  axis  at  the  center  of  mass  was 


90 

calculated  via  Equation  2  . 


[2] 


where  lzz  is  the  moment  of  inertia  about  the  oscillating  axis  of  the  pendulum,  R  is  the  distance 
from  each  wire  connection  to  the  center  of  the  axis  of  rotation,  m  is  the  mass  of  the  object,  g  is 
acceleration  due  to  gravity,  L  is  the  length  of  the  wires,  and  x  is  the  period  of  oscillation.  When 
adding  an  object  of  unknown  inertia,  the  mass  and  inertia  from  Equation  2  can  be  split  into 
components  of  the  frame  and  object: 

m  =  [3] 

*B3  ~  ^ BSS  ^  ^ OBB  [4] 


where  subscripts  P  and  O  represent  the  platform  and  object,  respectively.  Therefore,  the  moment 
of  inertia  of  the  object  can  be  calculated  using  Equation  520: 


System  Validation 

To  validate  the  system  set-up,  a  rectangular  aluminum  block  with  a  known  mass  (2.45 
kg)  and  known  moment  of  inertia  about  its  horizontal  and  vertical  axes  was  tested  in  both 
positions  using  the  same  protocol  that  was  used  to  test  the  RSP’s.  Data  was  collected,  and 
moments  of  inertia  were  calculated  and  compared  to  the  object’s  known  moments  of  inertia. 
System  error  about  each  axis  of  rotation  was  considered  as  the  difference  between  the  measured 
and  known  inertial  value  about  that  axis. 

Malalignment  between  the  object’s  and  the  pendulum’s  center  of  mass  can  introduce 
error  in  the  moment  of  inertia  estimations  as  well.  Malalignment  can  either  be  of  a  linear  nature, 
i.e.  if  the  object’s  center  of  mass  is  not  aligned  vertically  over  the  platform’s  center  of  mass,  or  it 
can  be  of  rotational  nature,  i.e.  the  object’s  principal  axis  is  not  vertical.  Great  care  was  taken  to 
align  the  prostheses  centers  of  mass  and  principal  axes  with  the  center  of  mass  and  axis  of 
rotation  (vertical)  of  the  triangular  platform;  however,  it  is  possible  that  these  centers  of  mass 
and  axes  were  not  perfectly  aligned.  Linear  malalignments  will  cause  a  shift  in  the  object  + 
platform  system’s  resultant  center  of  mass  and  the  pendulum  will  no  longer  oscillate  about  the 
intended  axis  of  rotation,  thus  imparting  error  in  the  period  and  moment  of  inertia  calculations. 
Rotational  malalignments  will  estimate  the  moment  of  inertia  about  an  axis  different  than  the 
axis  of  interest.  In  trifilar  pendulums,  it  is  often  impossible  to  ensure  exact  alignment  of  the 
centers  of  mass  and  axes  of  interest  of  the  system,  especially  when  using  non-uniformly  shaped 
objects  such  as  prostheses.  In  order  to  account  for  these  errors,  moment  of  inertia  was  measured 
at  varying  degrees  of  both  linear  and  angular  malalignment.  To  test  the  effects  of  linear 
malalignment,  the  aluminum  block,  in  both  the  horizontal  and  vertical  positions,  was  placed  1 
cm  to  10  cm  from  the  pendulum’s  axis  of  rotation  at  1cm  increments.  Moment  of  inertia  was 


calculated  following  the  same  protocol  described  earlier  and  compared  to  the  known  moment  of 
inertia  of  the  block  about  the  axis  of  interest.  Mathematically,  moment  of  inertia  error  due  to 
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linear  malalignment  can  be  expressed  by 

*  — [6] 
which  is  influenced  both  by  the  increased  inertia  of  the  system  (first  term)  and  the  change  in 
weight  distribution  and  center  of  rotation  due  to  the  new  resultant  center  of  mass  position 
(second  term,  in  parentheses)  caused  by  the  malalignment  of  the  plate  and  object’s  centers  of 


mass. 


*****  INSERT  FIGURE  5  ABOUT  HERE  ***** 


Rotational  malalignment  can  also  affect  the  accuracy  of  inertial  calculations  when 
oscillating  the  prosthesis  about  its  principal  axes.  Each  RSP’s  curved  design  prevented  it  from 
naturally  balancing  with  its  principal  axes  aligned  with  the  pendulum’s  axis  of  rotation.  To  test 
the  effects  of  rotational  malalignment,  the  middle  category  of  each  type  of  RSP  was  re-tested 
with  its  x-  and  z-axes  tilted  ±5°  (Figure  5).  Moment  of  inertia  was  calculated  following  the  same 
protocol  and  compared  to  the  moment  of  inertia  of  the  prosthesis  with  its  principal  axes  aligned 
to  vertical. 


Results 

Each  prosthesis’  measured  and  estimated  center  of  mass  positions,  including  the 
predicted  r  and  6  values  are  presented  in  Table  2.  Primary  axis  moments  of  inertia  are  presented 
in  Table  3.  For  all  prosthesis  designs,  the  y-axis  moment  of  inertia  corresponding  to  the 
anatomical  flexion/extension  axis  resulted  in  the  largest  moments  of  inertia  whereas  the  z-axis 
(anatomical  intemal/extemal  rotation)  had  the  smallest  moments  of  inertia. 


*****  INSERT  TABLE  2  ABOUT  HERE  ***** 


*****  INSERT  TABLE  3  ABOUT  HERE  ***** 

System  validation  with  the  aluminum  block  of  known  moment  of  inertia  showed  that  the 
trifdar  pendulum  was  accurate.  The  error  in  moment  of  inertia  was  -6.21  xlO"5  kg  m2  for  the 
horizontal  and  -2.65  x  10'6  kg-m2  for  the  vertical  position  of  the  aluminum  block,  which 
represent  a  1%  and  0.1%  error  in  the  results  respectively. 

Figure  6  shows  errors  in  moment  of  inertia  due  to  linear  malalignment.  Errors  are 
expressed  both  mathematically  (errors  in  moment  of  inertia  predicted  by  Equation  6, 
“calculated”)  and  as  directly  measured  using  the  trifilar  pendulum  (“measured”).  The  calculated 
and  measured  error  values  show  good  agreement  in  linear  malalignments  less  than  or  equal  to 
five  centimeters.  Effects  of  angular  malalignment  about  the  x-  and  z-axes  are  presented  in  Table 
4. 

*****  INSERT  TABLE  4  ABOUT  HERE  ***** 

*****  INSERT  FIGURE  6  ABOUT  HERE  ***** 


Discussion 

Mass,  center  of  mass  positions,  and  moments  of  inertia  were  estimated  from  three 
different  stiffness  categories  for  each  of  four  different  running-specific  prosthesis  designs. 
Variations  in  each  of  these  parameters  were  identified  due  to  different  prosthesis  designs  and 
stiffness  categories. 

The  center  of  mass  positions  for  each  prosthesis  were  calculated  using  a  reaction  board, 
and  a  predictive  equation  was  developed  to  estimate  these  positions  using  the  relationship 
between  the  most  proximal  point,  most  distal  point,  and  the  center  of  mass  position  of  a 


particular  prosthetic  design.  Using  the  ratio  between  the  head-COM  and  head-toe  vectors  along 
with  the  angle  between  these  vectors,  as  illustrated  in  Table  3,  the  equation  exactly  predicted  the 
centers  of  mass  when  the  category-specific  ratio  and  angles  were  used.  When  using  the  average 
ratio  and  angle  for  each  prosthesis,  a  majority  of  the  centers  of  mass  were  predicted  to  within  less 
than  one  centimeter  of  error  along  the  x-  and  z-axes.  Several  of  the  predictions,  however, 
resulted  in  greater  than  one  centimeter  of  error  with  a  maximum  error  of  2.8  centimeters.  The 
larger  errors  were  a  result  of  a  greater  range  of  ratio  and  angle  values  across  the  stiffness 
categories  of  a  particular  prosthetic  design.  Using  different  body  segment  parameter  models 
(e.g.  cadaveric  vs  body  scanning-based  regression  equations)  to  predict  the  intact  foot  center  of 
mass  position  are  shown  to  vary  by  greater  than  2  cm  in  the  predicted  positions  .  This  suggests 
that  center  of  mass  predictions  within  the  range  reported  in  this  study  are  reasonable;  however,  it 
is  recommended  that  direct  measurements  of  the  center  of  mass  position  within  prostheses  be 
used  when  possible  to  ensure  the  most  accurate  data  and  to  reduce  the  possible  errors  these  data 
will  induce  in  joint  kinetic  estimations  from  inverse  dynamics  equations. 

Intact  limb  inertial  property  measurements  can  vary  substantially  depending  on  the 
method  used  to  measure  these  parameters.  Rao  et  al.  (2006)  compared  segment  inertial  property 
estimations  using  six  methodologies  including  one  geometric  model,  two  cadaveric-based 
models,  and  three  mass  scanning  models  from  live  subjects.21  They  identified  significant 
differences  between  methods  in  all  inertial  properties  for  each  of  the  foot,  shank,  and  thigh 
segments.  Estimated  moments  of  inertia  for  intact  limbs  are  reported  to  differ  by  between 
0.0025-0.003 1  kg-m2  for  the  foot,  0.009-0.020  kg-m2  for  the  shank,  and  0.026-0.085  kg-m2  for 
the  thigh  depending  on  the  calculation  method22, 23 .  Moment  of  inertia  estimation  errors  induced 
in  the  current  study,  by  purposely  misaligning  the  aluminum  block  center  of  mass  with  the 
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pendulum  center  of  mass  during  oscillation  trials,  resulted  in  errors  up  to  0.008  kg-m  with  a 
10cm  malalignment  (Figure  6).  Malalignments  of  5cm  or  less  resulted  in  errors  less  than  0.002 
kg-m  ,  which  is  lower  than  differences  in  intact  foot  moments  of  inertia  due  to  measurement 
technique  differences.  This  indicates  that  using  a  trifilar  pendulum  to  estimate  moments  of 
inertia  of  prosthetic  components  will  yield  errors  less  than  those  currently  accepted  in  the 
literature  for  intact  limbs,  as  long  as  the  center  of  mass  of  the  prosthetic  device  is  aligned  within 
5cm  of  the  pendulum’s  center  of  mass.  With  careful  attention,  any  researcher  should  be  able  to 
reasonably  align  the  centers  of  mass  of  a  prosthesis  and  pendulum  with  errors  much  smaller  than 
5  cm,  which  would  yield  even  smaller  errors  in  moment  of  inertia. 

The  center  of  mass  and  moment  of  inertia  measurements  are  limited  to  the  prostheses  in 
their  uncompressed  form  (with  no  load  as  in  swing  phase  of  running).  The  inertial  parameters 
will  change  when  the  prosthesis  is  compressed  (for  example,  during  stance  phase  of  running). 
However,  since  the  loads  required  to  compress  the  prosthesis  are  very  large  relative  to  the 
changes  in  inertial  properties  during  loading,  it  is  likely  that  these  inertial  changes  would  have  a 
negligible  effect  on  the  resultant  inverse  dynamics  estimations  of  joint  kinetic  values.  Rather  the 
loads  (i.e.  ground  reaction  forces)  would  dominate  the  inverse  dynamics  predictions.  Additional 
studies  are  needed  to  determine  these  effects  and  to  discriminate  between  the  effects  of  inertial 
changes  and  the  loads  required  to  produce  those  changes. 

Conclusions 

The  inertial  properties  of  four  commonly  prescribed  running-specific  prostheses  were 
measured  using  a  reaction  board  and  a  trifilar  pendulum.  The  inertial  parameters  were  shown  to 
vary  slightly  between  stiffness  categories  within  a  prosthetic  design,  and  they  varied  more 
substantially  between  different  prosthetic  designs.  A  predictive  equation  was  presented  to 


estimate  the  center  of  mass  of  a  prosthesis  when  direct  measurements  are  not  possible.  These 
data  may  be  used  for  predicting  inertial  parameters  of  similar  prostheses.  The  predictive  equation 
and  pendulum  method  presented  measured  inertial  properties  with  errors  equal  to  or  less  than 
those  found  in  commonly  used  predictive  methods  for  intact  limb  inertial  parameters. 
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List  of  Table. 


Table  1.  Stiffness  categories  used  for  each  prosthesis  during  testing  and  the 
manufacturer  recommended  body  mass  range  associated  with  each  category. 

Prosthesis  Stiffness  Category  (body  mass  range) 


Flex-Run 

3  (53-59  kg) 

5  (69-77  kg) 

7  (89-100  kg) 

Nitro 

3  (60-68  kg) 

6  (89-100  kg) 

7  (101-116  kg) 

Cheetah 

3  (60-68  kg) 

5  (78-88  kg) 

7  (101-116  kg) 

1E90 

140  lb  (63.6  kg) 

185  lb  (84.1  kg) 

235  lb  (106.8  kg) 

Table  2.  Center  of  mass  (COM)  positions,  in  m,  along  the  principal  axes  measured  with  a  reaction 
board  in  the  sagittal  plane  (x-z  plane)  relative  to  the  “head”  position  (most  proximal  point  on  the 
prosthesis)  compared  to  COM  estimated  using  equation  1  in  the  text.  The  y-position  of  the  COM  is 
aligned  with  the  midline  of  the  prosthesis  and  thus  has  a  zero  value.  The  r  and  0  values  specific  to 
each  prosthesis  exactly  predicted  the  measured  COM.  The  average  r  and  0  values  measured  across 
stiffness  categories  for  a  particular  prosthesis  design  were  used  as  the  input  variables  to  predict  the 
estimated  COM  positions. 


Prosthesis  Cat 

.  „  .  „  Measured 

r  Avg  r  0  Avg  0  CQM 

Estimated 

COM 

Error 

X  z 

x  z 

X 

3 

0.398 

0.951 

-0.065 

-0.072 

-0.069 

-0.070 

0.004 

-0.002 

Flex-Run 

5 

0.418 

0.404 

1.044 

0.999 

-0.065 

-0.079 

-0.059 

-0.079 

-0.006 

0.001 

7 

0.395 

1.003 

-0.058 

-0.078 

-0.059 

-0.080 

0.001 

0.003 

3 

0.366 

1.128 

-0.053 

-0.069 

-0.062 

-0.054 

0.009 

-0.015 

Nitro 

6 

0.359 

0.362 

1.091 

1.185 

-0.051 

-0.069 

-0.063 

-0.053 

0.012 

-0.016 

7 

0.362 

1.336 

-0.058 

-0.063 

-0.054 

-0.061 

-0.004 

-0.002 

3 

0.504 

0.328 

0.025 

-0.264 

0.018 

-0.292 

0.007 

0.028 

Cheetah 

5 

0.589 

0.556 

0.377 

0.360 

0.019 

-0.304 

0.022 

-0.286 

-0.003 

-0.018 

7 

0.574 

0.374 

0.021 

-0.300 

0.024 

-0.290 

-0.003 

-0.010 

1401b 

0.597 

0.441 

0.017 

-0.307 

0.019 

-0.286 

-0.002 

-0.021 

1E90 

1851b 

0.538 

0.562 

0.452 

0.429 

0.029 

-0.278 

0.037 

-0.287 

-0.008 

0.009 

2351b 

0.537 

0.392 

0.032 

-0.278 

0.023 

-0.289 

0.009 

0.011 

Cat  =  stiffness  category 


r  =  ratio  of  measured  prosthesis  head-COM  to  head-toe  vector  magnitudes 
0  =  angle,  in  radians,  between  measured  head-COM  and  head-toe  vectors 
Error  =  difference  between  measured  and  estimated  COM  positions,  in  m 


Table  3.  Mass  and  moments  of  inertia  calculated  about  each  prosthesis’  measured 
primary  axis.  Category  represents  the  stiffness  category  of  the  prosthesis  according 
to  the  manufacturer. 

Moment  of  Inertia  (kg*m2) 


Prosthesis  Type  Category  Mass  (kg)  x  -  axis  y  -  axis  z  -  axis 


l  Ossur 

J  Cheetah 

3 

5 

7 

0.492 

0.511 

0.539 

0.0123 

0.0127 

0.0136 

0.0139 

0.0143 

0.0152 

0.0021 

0.0023 

0.0022 

3 

0.416 

0.0037 

0.0047 

0.0014 

1  Ossur 

5 

0.437 

0.0037 

0.0051 

0.0017 

y  t  lex-Kun 

7 

0.466 

0.0040 

0.0054 

0.0017 

3 

0.307 

0.0021 

0.0029 

0.0010 

A  Freedom 

J  Innovations 

6 

0.349 

0.0024 

0.0033 

0.0012 

Nitro 

7 

0.366 

0.0026 

0.0036 

0.0012 

1401b 

0.543 

0.0116 

0.0130 

0.0027 

Ottnhnrk 

1E90 

1851b 

0.605 

0.0131 

0.0152 

0.0035 

2351b 

0.677 

0.0144 

0.0168 

0.0042 

Table  4.  Error  values  of  angular  malalignment  of  ±5°  orientation.  Values  represent  the 
difference  between  the  prosthesis’  measured  moment  of  inertia  in  neutral  alignment  and 
±5°  rotation.  Data  are  presented  for  the  middle  stiffness  category  for  each  prosthesis. 


Prosthesis 


x-axis 


z-axis 


Flex-Run 

Nitro 

Cheetah 


-0.000009 

0.000654 

-0.000303 

0.000299 


0.000124 

0.000707 

0.000765 

0.001651 


0.000184 

0.000855 

0.000461 

0.000755 


0.000138 

0.000505 

-0.000502 

-0.000750 
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Figure  1.  Relationship  between  the  center  of  mass  (COM),  Head,  and  Toe  of  a 
running-specific  prosthesis.  T  is  the  Head-Toe  vector,  CM  is  the  Head-COM  vector, 
and  0  is  the  angle  between  these  vectors.  Equation  1  in  the  text  may  be  used  to 
estimate  the  COM  position  based  on  a  known  0  value  and  ratio  between  CM  and  T  for 
a  particular  running-specific  prosthesis  design. 


Figure  2.  Custom-built  trifilar  pendulum.  A  triangular  platform  is  suspended  from  a  large  frame 
by  three  equidistant  wires  that  allow  rotation  about  the  platform’s  center  of  mass.  The  moments 
of  inertia  of  a  prosthesis  may  be  calculated  directly  from  the  period  of  oscillation  of  the 
pendulum,  measured  when  one  comer  of  the  platform  passes  a  laser  sensor  (not  shown). 


a.  Freedom  b.  Ossur  Flex-Run  c.  Ossur  Cheetah  d.  Otto  Bock 
Innovations  Nitro  Sprinter  1E90 

Figure  3.  The  local  coordinate  systems  originating  at  the  centers  of  mass  for 
each  prosthesis.  The  y-axis  is  orthogonal  to  the  sagittal  plane,  with  the  positive 
direction  pointing  at  the  reader.  Details  on  defining  all  three  axes  are  included 
in  text. 


Axis  of  Rotation 


Figure  4.  Experimental  setup.  Primary  axes  of  rotation  (z-axis  shown  here)  of  the  prosthesis  are 
aligned  with  the  trifilar  pendulum’s  axis  of  rotation.  The  centers  of  mass  of  the  two  systems  are 
aligned. 


+5’ 


Figure  5.  Inducing  a  rotational  malalignment  of  the  z-axis  for  +/-  5°  tilt  of  the  intended 
prosthesis  axis  of  rotation. 
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♦  Horizontal:  Measured 
■  Vertical:  Measured 
A  Horizontal:  Calculated 
X  Vertical:  Calculated 


Figure  6.  Offset  error  as  measured  from  the  linear  malalignment  testing  with  the  aluminum 
block  and  through  mathematical  calculations.  Measured:  The  block  was  tested  in  both  the 
horizontal  and  vertical  position  and  displaced  linearly  from  the  intended  axis  of  rotation  in  1  cm 
increments.  Error  of  MOI  is  the  difference  between  the  known  MOI  of  the  block  when  aligned 
with  the  axis  of  rotation  and  the  measured  MOI  of  the  offset  block. 


